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Abstract. We consider the effects of a cosmological constant on the dynamics of constant angular momentum discs orbiting 
Schwarzschild-de Sitter black holes. The motivation behind this study is to investigate whether the presence of a radial force 
contrasting the black hole's gravitational attraction can influence the occurrence of the runaway instability, a robust feature of 
the dynamics of constant angular momentum tori in Schwarzschild and Kerr spacetimes. In addition to the inner cusp near the 
black hole horizon through which matter can accrete onto the black hole, in fact, a positive cosmological constant introduces 
also an outer cusp through which matter can leave the torus without accreting onto the black hole. To assess the impact of 
this outflow on the development of the instability we have performed time-dependent and axisymmetric hydrodynamical sim- 
ulations of equilibrium initial configurations in a sequence of background spacetimes of Schwarzschild-de Sitter black holes 
with increasing masses. The simulations have been performed with an unrealistic value for the cosmological constant which, 
however, yields sufficiently small discs to be resolved accurately on numerical grids and thus provides a first qualitative picture 
of the dynamics. The calculations, carried out for a wide range of initial conditions, show that the mass-loss from the outer cusp 
can have a considerable impact on the instability, with the latter being rapidly suppressed if the outflow is large enough. 
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1. Introduction 

Relativistic accretion tori orbiting around stellar-mass black 
holes have been the subject of renewed interest over the last 
few years in connection with the different astrophysical scenar- 
■ ios where these objects are expected to form, such as the core 
collapse of a massive star leading to a "failed" supernova explo- 
sion (a collapsar), or in the catastrophic merger of two (unequal 
mass) neutron stars in a close binary system. However, thick ac- 
cretion discs are probably present at much larger scales as well, 
suiTounding quasars and other active galactic nuclei, and feed- 
ing their central supermassive black holes. One of the major 
issues about such systems concerns their dynamical stability. 
This has important implications on the most favoured current 
models for the central engines of 7-ray bursts, either coUap- 
sars or binary neutron star mergers, for long and short bursts, 
respectively (see, e.g. Meszaros 2002 for a recent review). 

Discs around black holes may suffer from a number 
of instabilities produced either by axisymmetric or by non- 
axisymmetric perturbations and further triggered by the pres- 
ence of magnetic fields. A type of instability that has been stud- 
ied in a number of works and that could take place when the 
discs are geometrically thick and axisymmetric is the so-called 
runaway instability (see Font & Daigne, 2002a; Zanotti et al. 
2003 and references therein). To appreciate the mechanism 



leading to the development of this instability, consider an in- 
viscid fluid torus with a vertical structure and internal pressure 
gradients orbiting around a black hole (either Schwarzschild 
or KeiT). If the fluid is non self-gravitating, it will be con- 
tained within isopotential surfaces which generically possess 
a cusp on the equatorial plane (Fishbone & Moncrief, 1976; 
Kozlowski et al. 1978; Abramowicz et al. 1978). As a result, 
material from the disc can accrete onto the black hole through 
the cusp as the result of small deviations from hydrostatic equi- 
librium. 

Any amount of matter lost by the disc and captured by 
the black hole will increase its mass (and angular momen- 
tum), resulting in a modification of the equipotential surfaces 
which may cause the cusp to move deeper inside the torus more 
rapidly than the inner edge of the torus. When this happens, ad- 
ditional disc material will be allowed to fall into the black hole 
in an increasingly accelerated manner leading to the runaway 
instability. 

Although this instability was first studied in the '80s 
(Abramowicz et al. 1983; Wilson, 1984), time-dependent hy- 
drodynamical simulations have been performed only recently, 
either with SPH techniques and pseudo-Newtonian poten- 
tials (Masuda & Eriguchi 1997; Masuda, Nishida & Eriguchi 
1998), or with high-resolution shock-capturing (HRSC here- 
after) techniques in general relativity (Font & Daigne, 2002a 
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Zanotti et al. 2003). These investigations have shown that, 
under the (idealised) assumption of constant specific angular 
momentum distributions, relativistic tori around Schwarzschild 
and Kerr black holes are generically unstable to the runaway in- 
stability, if non self-gravitating. The inclusion of more generic 
initial conditions, however, can disfavour the occurrence of the 
instability. Recently, Font & Daigne (2002b) (see also Daigne 
& Font, 2003) have shown through numerical simulations that 
the runaway instability is suppressed when a non-constant dis- 
tribution of the angular momentum is assumed for the torus 
(increasing as a power-law of the radius), a result which is in 
agreement with studies based on a recent perturbative analysis 
(Rezzolla et al. 2003a; 2003b). While a similar stabilizing ef- 
fect has been shown to be provided by the black hole if this 
is rotating (Wilson, 1984; Abramowicz et al. 1998), Masuda & 
Eriguchi (1997) were able to show that the inclusion of the self- 
gravity of the torus effectively favours the instability. Clearly, 
a final conclusion on the occurrence of this instability has not 
been reached yet and will have to wait for fully general rela- 
tivistic simulations. However, the increasingly realistic inves- 
tigations performed recently have addressed several important 
aspects and the prospects are that we may be close to reaching 
a detailed description of the dynamics of the instability. 

A further physical process acting against the instability and 
which has not been investigated so far, is provided by the exis- 
tence of a repulsive force pointing in the direction opposite to 
the black hole's gravitational attraction. Such a force could dis- 
turb and even balance the standard outflow of mass through the 
inner cusp, thus potentially suppressing the runaway instability. 
As suggested recently by Stuchlik et al. (2000), such conditions 
could arise naturally in a black hole spacetime with a positive 
cosmological constant, i.e. in a Schwarzschild-de Sitter space- 
time. In such a spacetime, in fact, a second cusp appears in the 
outer parts of the equilibrium tori, near the so-called "static ra- 
dius". Assuming a value for the relict cosmological constant 
A ^ 10~^^cm^^ as deduced from recent cosmological ob- 
servations of the vacuum energy density (Krauss 1998) and 
compatible with a sample of observational estimates provided 
by the analysis of a large number of high redshift supernovae 
(Perlmutter et al. 1999; Riess et al. 1998), Stuchlik et al. (2000) 
find that the location of this outer cusp for the largest stationary 
discs which can be built in a Schwarzschild-de Sitter space- 
time is at about 50 — 100 kpc for supermassive black holes 
with masses in the range ^ 1G^A/q — IG^A/q. As for the in- 
ner one, a slight violation of the hydrostatic equilibrium at the 
outer cusp would induce a mass outflow from the disc and away 
from the black hole, which could affect the overall dynamics of 
the torus. 

However, this is not the only way in which a cosmolog- 
ical constant could modify the dynamics of a disc orbiting 
around a Schwarzschild-de Sitter black hole. As argued by 
StuchKk et al. (2000), in fact, a cosmological constant could 
produce a sensible modification in the accretion processes onto 
primordial black holes during the very early stages of expan- 
sion of the Universe, when phase transitions could take place, 
and the effective cosmological constant can have values in 
many orders exceeding its present value (Kolb & Turner 1990). 
Furthermore, a positive cosmological constant could also re- 



sult into strong collimation effects on jets escaping along the 
rotation axis of the central black hole (Stuchlfk et al. 2000). 

The aim of this paper is to investigate one of these in- 
triguing possibilities through numerical simulations. More pre- 
cisely, we present a comprehensive study of the nonlinear hy- 
drodynamics of constant angular momentum relativistic tori 
evolving in a sequence of background Schwarzschild-de Sitter 
spacetimes with increasing black holes masses. Our study clar- 
ifies the dynamical impact of a mass outflow on the occurrence 
of the runaway instability in such relativistic tori. 

We note that our setup will not be an astrophysically realis- 
tic one. There are two important reasons for this. Firstly, given 
the present estimates for the value of the cosmological constant 
and for the masses of the supermassive black holes believed 
to exist in the centre of galaxies and active galactic nuclei, we 
are still lacking sufficient observational evidence that stationary 
thick accretion discs exist on scales of about 100 kpc. Secondly, 
even assuming that such objects are present within large galax- 
ies, numerical calculations would have to face the present com- 
putational limitations which make it extremely hard to simu- 
late accurately accretion discs with very low rest-mass densi- 
ties and over such length scales. As a result, we will adopt a 
value for the cosmological constant that is unrealistically high. 
This yields discs with radial extents that are sufficiently small 
to be evolved numerically with satisfactory accuracy and pro- 
vides a first qualitative description of the role that a cosmolog- 
ical constant could play on the dynamics of relativistic tori. In 
addition to this, the calculations reported here also offer a way 
of assessing how the self-gravity of the torus, which is basically 
contrasting the gravitational attraction of the black hole, could 
modify the overall inertial balance. This will provide a use- 
ful insight when fully relativistic calculations solving for the 
Einstein equations coupled to a self-gravitating matter source 
will be performed. 

The organization of the paper is as follows. In Sect. |2l 
we briefly review the main properties of relativistic tori in a 
Schwarzschild-de Sitter spacetime. Next, in Sect.|3]we present 
the hydrodynamics equations and the numerical methods im- 
plemented in our axisymmetric evolution code. The material 
presented in this Section is rather limited, since the details 
have previously been reported in a number of papers. The 
last part of this Section is devoted to a discussion of the ini- 
tial data we use for the simulations. The numerical results 
are then described in Sect. |3 and, finally. Sect. |5] contains our 
conclusions. Throughout the paper we use a space-like signa- 
ture (— , +, +, +) and a system of geometrized units in which 
c = G = 1 . The unit of length is chosen to be the gravitational 
radius of the black hole, = GM / (?, where M is the mass of 
the black hole. Greek indices run from to 3 and Latin indices 
from 1 to 3. 

2. Stationary configurations in a 
Scliwarzscliiid-de Sitter spacetime 

Building on a wide literature discussing equilibrium con- 
figurations of perfect fluid relativistic tori orbiting around 
Schwarzschild or Kerr black holes, Stuchlfk et al. (2000) have 
recently extended these results to the case of a Schwarzschild- 
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de Sitter black hole. In spherical coordinates {t, r, 0, (j)) the line 
element of this spacetime reads 



1 



2M 



2M 

r 

„2 



y 



M2 



dr^ +r^{de^ +sin^ed(l)^) , (1) 



where M is the mass of the black hole and the cosmological 
constant A is incorporated in the dimensionless parameter y 
defined as 

y = ^AAf 2 . (2) 

This parameter has to be smaller than a critical value y < 
1 /27 in order to produce static regions of the spacetime where 
equilibrium configurations can be found. We note that a nega- 
tive cosmological constant, corresponding to a Schwarzschild- 
anti de Sitter black hole, does not introduce new qualitative fea- 
tures in the development of the runaway instability when com- 
pared to a Schwarzschild spacetime and will not be considered 
here. In what follows we briefly recall the main features of sta- 
tionary configurations in a Schwarzschild-de Sitter spacetime, 
referring the reader to the work of Stuchlfk et al. (2000) for an 
exhaustive discussion. 

The perfect fluid with four-velocity it^ is described by the 
usual stress-energy tensor 



(3) 



where g^'^ are the coefficients of the metric O in which, how- 
ever, the black hole mass M could be a function of time to 
account for the mass accreted onto the black hole [cf. the dis- 
cussion of Eq. i23\ below]. The fluid variables e, p, p, and 
h = {e + p)/p aie the proper energy density, the isotropic 
pressure, the rest mass density, and the specific enthalpy, re- 
spectively. An equation of state (EOS) of polytropic type, p = 
np^ — pe{j — 1), completes the thermodynamical description 
of the fluid. Here, k is the polytropic constant, 7 is the adiabatic 
index and e — e/p—lis the specific internal energy. As shown 
by Kozlowski et al. (1978) (see also Fishbone & Moncrief, 
1976), the pressure gradients can balance the gravitational and 
centrifugal forces, allowing for the existence of stationary con- 
figurations of matter in non-geodesic circular motion and con- 
tained within closed "constant pressure" equipotential surfaces. 
Under the conditions of hydrostatic equilibrium and of axisym- 
metry (i.e. dt = ~ 0) the relativistic Euler equations for 
a fluid with four-velocity u" = (w*, 0, 0, u'^) take the simple 
Bernouilli-type form 



^r,9, (4) 



where W — W{r, 9) = \n{ut) is the effective potential, £ = 
—u^/ut is the specific angular momentum, and fi = u'^/u* 
is the coordinate angular velocity as measured by an observer 
near the static radius, where the spacetime geometry is very 
close to a flat one (Stuchlik et al. 2000). Note that an explicit 
relation exists between the angular velocity and the specific an- 
gular momentum, which is given hy ft = —(-{gtt/ g<j}<t>)- 



Once M and A have been prescribed, the explicit expres- 
sion for the potential W{r,9) in the Schwarzschild-de Sitter 
spacetime is simply given by (Stuchlik et al. 2000) 



W(r,9) = iln 



(1 - 2M/r - yrVM^) 



sin^ I 



9 • 2 
sm 



6 -{I- 2M/r - yr'^/]VP)P 



(5) 



It is apparent from Eq. (0} that the simplest and indeed best 
studied configurations are obtained when the distribution of 
specific angular momentum I is prescribed to be constant. In 
this case, which is also the one considered here, W{r^ 0) shows 
three local extrema in the equatorial plane, one more than in the 
case of a Schwarzschild spacetime (cf. Fig.^. Ordering these 
points with increasing values of r, the first extremum corre- 
sponds to the position of the inner cusp, where the equipo- 
tential surface has a self-crossing point in the (r, 9) plane. The 
second extremum corresponds to the position of the "centre" of 
the torus, r,„ax, where the internal pressure of the torus has its 
maximum. The third extremum, finally, marks the position of 
the outer cusp, r^o, which is not present in the Schwarzschild 
spacetime and represents the distinctive contribution of the 
non-zero cosmological constant (see Fig.Qfor a schematic di- 
agram). Note that for a configuration with constant l, a local 
extremum of W is also a point of vanishing pressure gradients 
[cf. eq. (0}]. As a result, the fluid motion on a circular orbit of 
radius rd, r-max or r^o is a purely geodetic one, with the spe- 
cific angular momentum being given by the Keplerian one 
for a point-like particle at that radius. In a Schwarzschild-de 
Sitter spacetime this is given by 



_ r^M{l-yr^/M^) 
= {r - 2M - yr^ / IVPy 



(6) 



If the black hole is a Schwarzschild one, then y = and the 
admittable values for a constant distribution of specific angular 
momentum that will lead to a torus of finite size are 



(7) 



where ^,„s = 3-\/6/2 and ^mb = 4 represent the Keplerian 
specific angular momentum at the marginally stable and at 
the marginally bound orbits, respectively. In the case of a 
Schwarzschild-de Sitter black hole, however, the allowed pa- 
rameter space is more complex since in this case one needs 
to take into account also the admittable values for y. While a 
detailed discussion of this issue is provided by Stuchlik et al. 
(2000), we here remind that for y ^ 0, a toroidal configura- 
tion of finite size is obtained for values of the specific angular 
momentum that satisfy 



<£< 



(8) 



where now ^^ns ; and ^ms.o represent the values of the spe- 
cific angular momentum at the inner and outer marginally sta- 
ble radii, respectively, and correspond to the local minimum 
and maximum of the Keplerian specific angular momentum 
given by Eq. (|6j (cf. Fig. Hence, ^,„s,i and ^ms,o provide 
the minimum and maximum values of i{r, y) for which stable 
Keplerian circular orbits exist. Note that the very existence of 
these stable circular orbits depends on the values of y and, in 
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Fig. 1. Schematic diagram for the Keplerian specific angular momentum Ij^ in a Schwarzschild-de Sitter spacetime (left panel) 
and the corresponding effective potential W (right panel) once a constant value for i has been chosen (cf. Fig|2j. The figure 
reports the different radial locations that are relevant for our discussion: the inner and outer cusp points rci, Tco, the inner and 
outer radii for the torus rin,rout, the inner and outer marginally stable orbits rms.i, ^'ms.o, and the location of the maximum 
pressure in the torus r^ax (see text for details). Note that rci, fco and r^ax are determined once a value for the constant specific 
angular momentum has been chosen (this is shown with the long-dashed line in the left panel) and that the inner and outer radii 
need not coincide with the corresponding locations of the cusps but are set by the value chosen for the potential W-^ (this is shown 
with the long-dashed line in the right panel). Reported for comparison with a short-dashed line are and in a Schwarzschild 
spacetime; the radial coordinate is shown on a logarithmic scale. 



particular, these can be found only if y < ^ms = 12/15^ ^ 
2.37 X 10-^ < Ve- 
in other words, given a value for the dimensionless param- 
eter y < i/ins, circular orbits can exist for ^ms.i(y) < t < 
^ms.o(y)- Clearly, these orbits will be at radii smaller than the 
static radius Vg = at which the angular momentum (|6jl 

is zero and where the gravitational attraction of the black hole 
is exactly balanced by the cosmological repulsion'. 

Finally, it is worth reminding that for point-like particles, 
the decreasing part of £^ (i.e. for r < r^^^y, r > r^g □) cor- 
responds to unstable circular geodesies, while the increasing 
part (i.e. for rms.i < r < r^ns.o) corresponds to stable circu- 
lar geodesies. For a perfect fluid, however, internal pressure 
gradients allow the torus to occupy a region which would be 
geodesically unstable for point-like particles. A schematic dia- 
gram showing the behaviour of the Keplerian specific angular 
momentum and the corresponding effective potential W in 
a Schwarzschild-de Sitter spacetime, as well as the locations of 
the different radii that have been discussed so far is presented in 
Fig.n The Schwarzshild case is also displayed for comparison 
with a dashed line. 



' The Keplerian specific angular momentum is not well defined for 
r > rs, where its square is negative. 



Once a value for the cosmological parameter and for the 
angular momentum have been fixed, Eq. 0, supplemented by 
the polytropic EOS, can be integrated analytically for any r < 
rs, to yield the rest-mass density distribution inside the torus 

r7-i ^1/(7-1) 

p(r, 9) = I i— [exp(Win -W)-\]\ , (9) 

where Win = W^(nn,7r/2) and rjn is the inner edge of the 
torus. The latter is assumed to be a free parameter and is effec- 
tively controlled by the potential gap AWi = W{rin, 7r/2) — 
W{r,u7r/2). 

Hereafter, we will focus on tori built in a parameter space 
that is smaller than the one discussed so far In particular, we 
will consider tori with constant specific angular momentum in 
the range e„,sA < i < iph < 4is,o, where ^^i^ = r^/{r - 
2M — yr^ /AP) is the angular momentum of the unstable 
photon circular geodesic (Stuchlik et al. 2000). Furthermore, 
the hydrodynamical evolution of these tori will be followed in 
Schwarzschild-de Sitter spacetimes with dimensionless cosmo- 
logical constant < y < ye 1/118125 ~ 8.46 x IQ-*^, 
where y^ corresponds to the value of y for which the minimum 
of £ph is equal to the Keplerian angular momentum of the outer 
marginally stable orbit (StuchKk et al. 2000). 
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3. Hydrodynamic Equations and Numerical 
Approachi 

3. 1. Basic Equations 

Suitable initial data for the torus (see Sect. l3.3l for specific de- 
tails) are evolved forward in time using a two-dimensional (2D) 
general relativistic hydrodynamics code. The equations for a 
perfect fluid are implemented using the formulation developed 
by Banyuls et al. (1997) in which, adopting a 3 + 1 decomposed 
spacetime, the covariant conservation laws of the density cur- 
rent and of the stress-energy tensor are recast into a first-order, 
flux-conservative system of the type 

'dt ' 



dr 



de 



(10) 



where a = \/—gm is the lapse function of the Schwarzschild- 
de Sitter metric and U = (£>, S'r, S'e, S'^) is the state-vector 
of "conserved" variables. These represent the relativistic rest- 
mass density and the relativistic momenta in the three co- 
ordinate directions, respectively. The conserved variables are 
expressed in terms of the "primitive " variables p, Vi , and e 
through the relations D = pF and Sj = phT^vj, where the 
covariant components of the three-velocity Vi can be obtained 
from the contravariant components = jau^ through the 
coefficients of the spatial 3-metric 7^, i.e. Vi — ^ijvK The 
Lorentz factor F is assumed to be measured by a local static 
observer and is given by F = au* = (1 — u^)^^/^, with 
= ^ijV^vK Finally, the vectors F' and S appearing in 
Eq. ( I10> represent the fluxes and sources of the evolved quanti- 
ties, respectively, the latter being due entirely to the spacetime 
curvature. The explicit expressions of these terms can be easily 
generalized from the Schwarzschild case and are 

F'-(w) = V') > (H) 

F\^) = {Dv\Srv\Sev'' ^V.S^v^) , (12) 

5(w) = 52, 53,54) , (13) 
where 

Si = DiAv-" - Bv^) , 



(14) 



S2 



^-4t? ) (r + D)-acot05,i;'' + 



a 



{Sev' + S^v"^ - 2Srv'') , (15) 



SeiAv'' - Bv^) + ar^S-^t;'^ sii 



54 = S^{Av'' — 

and where 

M 2a 



A 



ry 



B = acotf 



(16) 



(17) 



(18) 



(19) 



It should be noted that since the dynamical evolution is as- 
sumed to be adiabatic, no evolution equation for the relativistic 
total energy density t = phV"^ — p ~ D is solved here. 



3.2. Numerical approacli 

The numerical code we use is the same employed by Zanotti et 
al. (2003) to study the dynamics of constant angular momen- 
tum relativistic discs around a Schwarzschild black hole. For 
the present investigation the code has been extended to account 
for the modifications introduced by the Schwarzschild-de Sitter 
geometry. The general relativistic hydrodynamics equations are 
solved by means of a HRSC scheme based on Marquina's flux 
formula (see e.g. Font, 2000 for a review of these methods in 
numerical general relativistic hydrodynamics). 

In order to cover optimally the large spatial extent of the 
equilibrium configurations and yet reach a satisfactory spa- 
tial resolution in the regions closer to the two cusps where 
the fluid motion needs to be calculated most accurately, we 
have introduced an important technical modification in the han- 
dling of the radial-coordinate grid. More precisely, we use a 
non-uniform radial grid with a logarithmic spacing, which is 
double- varied in the vicinities of the inner and outer cusps. 
The coordinate mapping used for this purpose is reminiscent 
of a tortoise coordinate mapping but it has been extended to a 
Schwarzschild-de Sitter metric as 

= ± / ( 1 - — - y— I dr , (20) 



M2 / 



where the ± sign distinguishes whether the mapping is made 
for increasing or decreasing values of the coordinate r, respec- 
tively. 

As a result of this mapping, a radial grid of Nr = 300 
zones allows to cover a spatial domain going from r^^j^^ =2.1 
to r^j^j, = 100 with a minimum radial spacing of the innermost 
part of the grid Ar = 10~^ and, correspondingly, a minimum 
radial spacing Ar = 10~^ for the outermost part of the radial 
grid. The two grids join smoothly at r = 48.6, where the res- 
olution is Ar = 2.72. The angular grid, on the other hand, is 
more straightforward to build and consists of Ng — 70 equally 
spaced zones extending from to tt (cf. Zanotti et al. 2003) 

As in Zanotti et al. (2003), a low density "atmosphere" is 
introduced in those parts of the numerical domain outside the 
torus. The initial atmosphere model chosen corresponds to the 
spherically symmetric accretion solution of non-interacting test 
fluid particles. The maximum density of the atmosphere is typ- 
ically 5 to 6 orders of magnitude smaller than the density at the 
centre of the torus. In all of the validating tests performed, the 
hydrodynamical evolution of the torus was found to be unaf- 
fected by the presence of this atmosphere, which is evolved as 
the bulk of the fluid. 

Finally, the mass outflows at the innermost and outermost 
radial points are computed respectively as 



-2tt 



and 



Wo(Vax) = 27r 



^gDv^'de 



-gDv^'de 



(21) 



(22) 



where g is the determinant of the metric and \/—g = sin 6. 

Note that the mass outflow given by Eq. J21l i corresponds 
effectively to the mass accretion rate onto the black hole and 
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Table 1. Main properties of the tori considered in the numerical calculations. From left to right the columns report: the name of 
the model, the specific angular momentum I (normalized to M), the polytropic constant k, the inner and outer cusps of the torus, 
and Tco, the radial position of the pressure maximum r^ax (^H radii are in units of the gravitational radius 7-g), the potential 
gaps AWi = Win — Wei and AWq = Win — W^o, where Win is the potential at the inner edge of the disc. The last column 
reports the orbital period at the centre of the torus, ^orb^ expressed in milliseconds. All of the models share the same value of the 
cosmological parameter y = 10~^, the same mass for the black hole, M — IOA/q, the same adiabatic index 7 = 4/3, and the 
same torus-to-hole mass ratio Mt/M = 0.2. 



Model 


I 


K (Cgs) 




^CO 


^max 


AWi 


AWo 


^orb (ms) 


Ai 


3.84 


8.970x10'" 


4.419 


94.866 


8.822 


0.010 


-0.010 


8.11 


A2 


3.84 


2.568x10^'^ 


4.419 


94.866 


8.822 


0.025 


0.005 


8.11 


Az 


3.84 


4.372x10^'^ 


4.419 


94.866 


8.822 


0.032 


0.012 


8.11 


Bi 


3.94 


2.295x10'" 


4.133 


94.564 


9.876 


0.004 


0.004 


9.61 


B2 


3.94 


3.775x10^" 


4.133 


94.564 


9.876 


0.010 


0.010 


9.61 


B2 


3.94 


6.740x10^" 


4.133 


94.564 


9.876 


0.020 


0.020 


9.61 


Ci 


4.00 


3.025x10''' 


4.000 


94.373 


10.489 


-0.007 


0.007 


10.51 


C2 


4.00 


7.120x10^'^ 


4.000 


94.373 


10.489 


0.007 


0.021 


10.51 


C2 


4.00 


1.125x10^'' 


4.000 


94.373 


10.489 


0.020 


0.034 


10.51 



is used to account for the instantaneous increase of the black 
hole mass at every time step. This, in turn, provides infor- 
mation about the changes in the background spacetime, fun- 
damental for the appearance of the runaway instability (Font 
& Daigne, 2002a; Zanotti et al. 2003). As mentioned in the 
Introduction, since we neglect the self-gravity of the torus, 
the hydrodynamics equations are solved in a sequence of 
background Schwarzschild-de Sitter spacetimes with increas- 
ing black hole masses. In practice, the spacetime evolution is 
achieved through a remapping of the metric functions at each 
time level of the type 

5^,(r,M",y)-^5^,(r,M"+\2/), (23) 

where = A/" + M mp(r^jj^) is the mass of the black 

hole at the new timelevel A detailed discussion on the 

validity of this approximation can be found in Font & Daigne 
(2002a) and in Zanotti et al. (2003). The prescription ( I23> is 
justified and can be regarded as a very good approximation 
when the variation of the black hole mass per unit time, m", 
is very small. This is certainly the case for the small disc-to- 
hole mass ratios Mt/M considered here (cf. Table Q. 

3.3. Initial Data 

As shown by Font & Daigne (2002a) and Zanotti et al. (2003), 
the runaway instability is a robust feature of constant angular 
momentum relativistic tori in Schwarzschild and Kerr space- 
times, if these are non self-gravitating. This result does not de- 
pend on the way the instability is triggered, i.e. by either ar- 
tificially expanding the torus over the potential barrier at the 
inner cusp (Font & Daigne, 2002a), or by introducing pertur- 
bations in an otherwise stable torus (Zanotti et al. 2003). As 
discussed above, however, the presence of an outer cusp in a 
Schwarzschild-de Sitter spacetime is likely to affect the robust- 
ness of this conclusion. In order to investigate to what extent 
an outflow of mass can interfere with the development of the 



runaway instability we have studied the behaviour of three dif- 
ferent classes of models, which we refer to as A, B and C. 
These models are distinguished on whether the effective po- 
tential at the inner cusp, Wei, is less than, equal to, or larger 
than the effective potential Wco at the outer cusp, respectively. 
Furthermore, for each of these classes of models we have con- 
sidered three different initial configurations, as shown in Fig.|2j 
with the potential at the inner edge of the torus, Wn, being dif- 
ferent from the potential barrier, Wei, by the adjustable amount 
AWi. 

Note that for all of the models of the class A, the hy- 
drostatic equilibrium is always violated at the inner cusp, i.e. 
Win > Wei, and a mass outflow will necessarily take place 
at the inner edge of the disc once the initial data is evolved. 
Furthermore, a mass loss will take place also at the outer cusp 
for models A2 and A3, which have Wn > Weo > Wei. For all 
of the models of the class B, on the other hand, the hydrostatic 
equilibrium is violated at both cusps and by the same amount, 
i.e. Win > Wei = Weo and, again, the mass outflows are regu- 
lated by the potential jump AW. Finally, for all of the models 
of the class C, the hydrostatic equilibrium is always violated at 
the outer cusp, i.e. Wn > Weo, and also at the inner cusp for 
models C2 and C3, for which Win > Wei > Wcq. All of these 
different initial conditions are illustrated in Fig.|2l whose dif- 
ferent panels show the effective potential curves and the values 
of Win for the three classes of models. 

Summarized in Table [2 are the potential jumps AW and 
A Wo, as well as the the most relevant parameters of the various 
initial models considered here. Note that we have used a poly- 
tropic index 7 = 4/3 and adjusted the polytropic constant so 
as to have a small torus-to-hole mass ratio Mt/M = 0.2, thus 
minimizing the error introduced by neglecting the self-gravity 
of the torus. 

A final comment will be made on the value adopted for the 
dimensionless cosmological constant y. As mentioned in the 
Introduction, the primary goal of our investigation is to pro- 
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Fig. 2. Radial profiles of the potential at the equatorial plane for models labelled A (left), B (centre) and C (right), respectively. 
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edges of the torus for each of the models considered. Models A, B and C differ among each other by the value of the angular 
momentum I which is reported in the lower right corners of the different panels. 



vide a first qualitative analysis of the dynamics of thick discs 
in spacetimes with a positive cosmological constant. Because 
of this, and in order to avoid shortcomings with observations or 
computational resources, we have performed numerical simu- 
lations with y — 10^^ for a black hole of mass M = IQMq. 
These parameters, although unrealistic, allow us to construct 
sufficiently compact discs and perform numerical simulations 
with a satisfactory spatial resolution. We expect that all of the 
quaUtative conclusions we draw from our simulations will con- 
tinue to hold when considering more realistic values for both 
the cosmological constant and the black hole masses. 

4. Results 

As mentioned above, the development of the runaway instabil- 
ity appears to be a robust feature of the dynamics of non self- 
gravitating tori orbiting around Schwarzschild or Kerr black 
holes with constant distributions of specific angular momen- 
tum. This conclusion has been reached after numerous simu- 
lations have been performed for a large range of torus-to-hole 
mass ratios Mt/M and under a number of different initial con- 
ditions (Font & Daigne, 2002a; Zanotti et al. 2003). In all of 
these simulations, the onset and full development of the insta- 
biUty was observed when the spacetime geometry was suitably 
modified to account for the black hole's mass-increase due to 
accretion. As a result of the instability, the torus is very rapidly 
accreted onto the black hole and this is most clearly signalled 
by the exponential growth in time of the rest-mass accretion 
rate at the innermost radial grid point. The growth-time for the 
instability is inversely proportional to the ratio M^jM and is 
comparable with the dynamical (i.e. orbital) timescale when 
MtjM - 1. 

In the case of a Schwarzschild-de Sitter black hole, how- 
ever, the rest-mass of the torus can change not only because 
of losses through the inner cusp leading to accretion onto the 



black hole, but also because of outflows from the outer cusp 
and away from the black hole. While both allowed, the impact 
that these two mass outflows could have on the dynamics of the 
torus is very different. The first one, in fact, induces changes in 
the black hole mass and could therefore lead to the runaway 
instability. The second one, on the other hand, does not pro- 
duce changes of the background spacetime and cannot there- 
fore produce an instability. Nevertheless, it can affect the hy- 
drodynamical evolution in a number of different ways: firstly, 
by reducing the amount of rest-mass in the torus available for 
accretion and, secondly, by producing significant alterations of 
the velocity field, especially in the outer regions of the torus. 

We have followed the hydrodynamical evolution of the 
models described in Tabled over a number of orbital periods 
sufficiently large to reveal the impact of a positive cosmolog- 
ical constant on the occurrence of the runaway instability. In 
Figs. |3] 13 and we show the two mass outflows rn\ and ttIq 
as a function of the orbital period <oib — STr/f^max at the cen- 
tre of the torus, and for the three classes of models listed in 
Tabled The three small insets shown in the panels for mj offer 
a view of the evolution of the rest-mass of the torus after this 
has been normalized to its initial value. The description of the 
dynamics of the tori is also completed with Figs.0]and|6l which 
show equally spaced isocontours of the logarithm of the rest- 
mass density and, superimposed, the velocity field for models 
A2 and i?2, respectively, at two different times during the evo- 
lution. 

As it is apparent from a rapid look at these figures, the run- 
away instability is no longer the only possible evolution of the 
system, whose dynamics is instead the result of the interplay 
between the inner and the outer mass outflows. The occurrence 
of the runaway instability is clearly visible in the left panel of 
Fig EJfor model A\ (solid line). Model A\, in fact, has initial 
conditions that resemble those encountered for a Schwarzschild 
black hole, with the outer radius of the torus located far from 
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Fig. 3. Time evolution of the inner (left panel) and of the outer (right panel) mass outflows for the models of class A. The data is 
shown in units of solar masses per second, while the time is expressed in units of the orbital period. Note that only models Ai 
(solid Une) and A2 (dotted line) are runaway unstable. The solid circles in the two panels indicate the time at which nii > rho 
for model A2 ■ 



the outer cusp (cf. Figs. ^|2jl. As a result, the right panel of 
Fig-El shows that the outer mass outflow is in this case very 
small (indeed slightly negative as a result of accretion onto the 
torus of the infalling atmosphere), while the mass accretion rate 
onto the black hole (left panel) grows exponentially and undis- 
turbed until the full development of the runaway instability at 
t ~ 24.2 toib- The dynamical evolution is different for mod- 
els A2 (dotted line) and ^3 (dashed line), where the competi- 
tion between the two outflows at the inner and outer edges of 
the disc is closer to a balance and the initial outer mass-loss 
is non-negligible. For model A2, in particular, this is clearly 
visible in the left panel of Fig. |3 which shows that at early 
times it ^ 1.2 iorb) the largest fluid velocities (w ~ 4 x 10^'^) 
are reached in the outer regions of the torus and are outwardly 
directed. However, the outer mass-loss eventually becomes in- 
sufficient to prevent the development of the runaway instability, 
which takes place after t ~ 33 forb- The corresponding veloc- 
ity field (with all vectors pointing towards the black hole) and 
the isocontours of the logarithm of the rest-mass density dis- 
played in the right panel of Fig.|3(at time t ^ 30.0 torb) show 
the important reduction in size undergone by the torus, which 
is about to disappear entirely inside the black hole after a few 
more orbital periods. This type of evolution does not take place 
for model A^, whose dynamics is completely dominated by the 
mass outflow through the outer cusp and for which the runaway 
instability does not develop (cf. Fig.|3}- 

To better interpret the dynamics behind these simulations it 
is useful to compare the amount of rest-mass in the torus after 
the first 10 orbital periods of the evolution for the three models 
of class A (see inset in the left panel of Fig. |3}. The residual 



rest-mass of the torus is 99%, 90%, and 55% of the initial one 
for the models Ai, A2, and A3, respectively. This may appear 
somewhat surprising given the fact that the mass outflow from 
the inner cusp is smaller in the case of model A^ than it is 
for model A2, despite A3 having a larger initial potential jump 
AWi. The explanation for this comes from looking at the right 
panel of Fig. |3l which shows that the mass outflow from the 
outer cusp is however larger for model A3 than it is for model 
A2. As a result, the torus is emptied more efficiently, and this 
happens mostly through the outer cusp. 

It is also very instructive to compare the mass outflows at 
the inner and outer edges of the disc for the three models of 
class A. Doing so leads to the important result that rhi > ttIq 
at all times for model Ai and that mj < rho at all times for 
model A3. In other words, the unstable and the stable models 
seem to differ from each other on whether the mass outflow 
from the inner edge is larger or smaller than the corresponding 
mass-loss from the outer edge. In the case of model A2, on the 
other hand, the two mass outflows are closer to a balance and 
TOi > mo only for t > 18 iorb, after which the runaway in- 
stability clearly develops; the time when this happens is shown 
with the filled circles in the two panels of Fig. |3l It appears, 
therefore, that an increasing potential gap at the outer edge of 
the torus/avoMri the outflow of mass from the outer cusp but it 
also disfavours (and to a larger extent) the mass outflow from 
the inner cusp, with the corresponding accretion onto the black 
hole being severely suppressed. As a result, the feed-back of 
the black hole spacetime onto the overall dynamics of the sys- 
tem is considerably reduced and the runaway instability sup- 
pressed. Stated differently, the development of the runaway in- 
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Fig. 4. Velocity field and equally spaced isocontours of the logarithm of the rest-mass density for model A2 at an early time 
(left panel) and at a later time (right panel); the times reported are in units of the orbital period. Initially the outer mass flux 
dominates the dynamics of the torus. However, the gravitational attraction of the black hole eventually overcomes the effect 
of the cosmological constant and the runaway instability takes place. This leads to the large inward-directed fluxes and to the 
disappearance of the torus inside the black hole in a few orbital periods. 
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Fig. 5. Same as Fig.|3]but for the models of class B. 



stability appears to be related to the efficiency of the mass-loss 
through the edges of the disc and, in particular, the instability is 
effectively "extinguished" whenever the mass outflow from the 
outer edge of the disc and away from the black hole is larger 
than flie mass outflow from the inner edge of the disc and onto 



the black hole. The condition rhi < rho can thus be used as a 
simple sufficient condition for the suppression of the runaway 
instability in a thick disc orbiting around a Schwarzschild-de 
Sitter black hole. 
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Fig. 6. Same as Fig.l^but for model i?2. The intense mass ouflow across the outer edge of the disc removes a large fraction of its 
mass, and suppresses the runaway instabiUty. The final disc reaches a quasi-steady state. 
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Fig. 7. Same as Fig.|3]but for the models of class C. 



The role played by a positive cosmological constant on the two panels of Fig.|5]show that the mass fluxes through the outer 
dynamics of the discs and described so far for the models of edges of the discs to infinity are always larger than the ones to- 
class A is present also for the models of class B and C, al- wards the black hole (i.e. mj < mo). As a result, the models of 
though with some slight differences. As discussed in Sect. 13.31 class B are all stable to the runaway instability. This is partic- 
and illustrated in Fig. |2] the models of class B are built with ularly apparent in models B2 and B^, for which the outer mass 
the outer and inner edges having the same effective potential, outflows are at least a couple of orders of magnitude larger 
One would therefore expect that this would yield to very similar than the corresponding mass outflows onto the black hole, and 
mass outflows at the two boundaries of the disc. However, the which become negligibly small (i.e. rhi < 10^*A/q/s) wefl 
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before 10 orbital periods. As a result, a large amount of the 
matter in those discs is not accreted onto the black hole, but es- 
capes to infinity. This is illustrated in the small inset of the left 
panel of Fig. |5] which shows that after about 20 orbital periods 
more than 60% of the torus rest-mass is lost for model B2 and 
more than 90% for model B3. 

Once the outflows die-off in the tori of class B, the remain- 
ing matter reaches a quasi-stable equilibrium, accreting onto 
the black hole on a timescale which is essentially controlled by 
the rate of mass-loss through the inner cusp. The importance 
of the mass outflow at the outer edges of the discs of class B 
is also apparent from Fig. |6l which shows the velocity field 
and isocontours of the logarithm of the rest-mass density of 
model B2 at an early and a later stage of the evolution. Note 
how the left panels of Fig. |4] and |6l have velocity fields that 
differ mostly in modulus but are equally oriented, while the 
right panels are substantially different with velocity fields that 
have opposite orientations leading to the disappearence of the 
torus into the black hole and to infinity, respectively. We also 
note that while the difference between the inner and outer mass 
outflows remains large also in the case of model Bi, the dy- 
namics is in this case much closer to an equilibrium, with the 
torus being still progressively emptied to infinity, but on a much 
larger timescale. No runaway instability was observed for this 
model over the time for which the calculations were carried out 

(t ~ 33 torb)- 

Finally, for models C (see Fig.Q, the dynamics of the discs 
is particularly simple and the final result is rather clear to inter- 
pret. In this case, in fact, all the discs are built with an effec- 
tive potential which is larger at the outer edge (cf. Fig.|2j and 
represent, therefore, initial conditions that are conceptually the 
opposite of those in models A. Because of the high potential 
barrier at the inner edge of the disc, the inner mass outflow is 
always rather minute and several orders of magnitude smaller 
than the corresponding mass outflow from the outer edge. As a 
result, the mass in the torus is lost very rapidly to infinity and 
very little is accreted onto the black hole. In particular, in the 
most dramatic case of model C3, the residual rest-mass in the 
torus is less than 20% after only 6 orbital periods. 

5. Conclusions 

We have investigated the effect of a positive cosmological con- 
stant on the dynamics of non self-gravitating thick accretion 
discs orbiting Schwarzschild-de Sitter black holes with con- 
stant distributions of specific angular momentum. The moti- 
vation behind this investigation has been that of assessing the 
role played by an effective repulsive force in the onset and de- 
velopment of the runaway instability, which represents a robust 
feature in the dynamics of constant angular momentum tori. In 
addition to the inner cusp near the black hole horizon, through 
which matter can accrete onto the black hole when small devi- 
ations from the hydrostatic equilibrium are present, thick discs 
in a Schwarzschild de-Sitter spacetime also possess an outer 
cusp through which matter can leave the torus without accret- 
ing onto the black hole. As a result of this mass-loss to infinity, 
the changes in the background metric (which are responsible 



for the development of the runaway instability) may be altered 
considerably and the instability thus suppressed. 

As a simple way to evaluate this effect we have consid- 
ered a sequence of Schwarzschild-de Sitter spacetimes differ- 
ing only in their total mass and have performed time-dependent 
general relativistic hydrodynamical simulations in these back- 
ground metrics of thick discs which are initially slightly out 
of hydrostatic equilibrium. In doing this we have adopted an 
unrealistically high value for the cosmological constant which 
however yields sufficiently small discs (extending up to about a 
few hundred gravitational radii) to be accurately resolved with 
fine enough axisymmetric numerical grids. 

We have performed a number of simulations involving 
initial configurations of constant specific angular momentum 
discs differing both for the relative amplitude of the peaks in 
the effective potential and for the potential jump at the inner 
and outer cusps. The results obtained indicate that the runaway 
instability is no longer the only possible evolution of these sys- 
tems but that their dynamics is rather the end-result of the in- 
terplay between the inner and the outer mass outflows. On the 
one hand, in fact, we have evolved initial models for which 
the cosmological constant has a weak influence; these mod- 
els have negligible mass outflows to infinity while maintaining 
large mass outflows onto the black hole, which then lead to the 
development of the runaway instability. On the other hand, we 
have evolved initial models which are significantly influenced 
by the cosmological constant; these models develop mass out- 
flows through the outer cusp which are much larger than those 
appearing at the inner cusp and, hence, do not develop the run- 
away instability. Placed somewhere between these two classes 
of initial configurations there exist initial models for which the 
mass outflows from the inner and outer cusps are more closely 
balanced. In these cases the runaway instability may or may 
not develop and we have noticed that a simple comparison be- 
tween the mass outflows can be used to deduce the fate of the 
accreting disc. More specifically, we have found that the con- 
dition rhi < trio provides a simple sufficient condition for the 
suppression of the runaway instability in a thick disc orbiting 
around a Schwarzschild-de Sitter black hole. 

In spite of the idealized setup used, the simulations per- 
formed here provide a first qualitative description of the com- 
plex nonlinear dynamics of thick discs in Schwarzschild-de 
Sitter spacetimes and we expect that most of the results ob- 
tained will continue to hold also when more realistic values 
for the cosmological constant are used. Aa a final comment we 
note that besides providing a qualitative description of the role 
that a cosmological constant could play on the dynamics of rel- 
ativistic tori, these calculations also offer a way of assessing, 
at least qualitatively, the inertial role that the self-gravity of the 
torus plays in the development of the runaway instability. This 
will be very useful when studying the dynamics of relativistic 
tori with numerical codes solving also the full Einstein equa- 
tions. 

Acknowledgements. It a pleasure to thank Z. Stuchlfk for use- 
ful comments. Financial support for this research has been pro- 
vided by the Italian MlUR, by the Spanish Ministerio de Ciencia 
y Tecnologi'a (grant AYA 2001-3490-C02-01) and and by the EU 



12 



Zanotti, RezzoUa and Font: Dynamics of thick discs around Schwarzschild-de Sitter black holes 



Network Programme (Research Training Network Contract HPRN- 
CT-2000-00137). LR also acknowledges the kind hospitality at the 
KIT? in Santa Barbara, where part of this research was carried 
out (NSF grant PHY99-07949). The computations were performed 
on the Beowulf Cluster for numerical relativity "AlbertlOO", at the 
University of Parma. 

References 

Abramowicz M. A., Jaroszynski M., Sikora M., 1978, A&A, 63, 221 
Abramowicz M. A., Calvani M., Nobili L., 1983, Nature, 302, 597 
Abramowicz M. A., Karas V., Lanza A., 1998, Ac&A, 331, 1143 
Banyuls, P., Font, J. A., Ibanez, J. M., Marti, J.M., Miralles, J. A., 

1997, ApJ, 476, 221 
Daigne, P., Pont, J. A., 2003, in preparation 
Fishbone L. G., Moncrief V., 1976, ApJ, 207, 962 
Font J. A., 2000, Living Reviews in Relativity, 3, 2 
Font J. A., Daigne F, 2002a, MNRAS, 334, 383 
Font J. A., Daigne R, 2002b, ApJL, 581, L23-L26 
Kozlowski M., JaroszyAski M., Abramowicz M. A., 1978, A&A, 63, 

209 

Kolb E.W., Turner M.S., 1990, The Early Universe. The Advanced 
Book Program, Addison- Wesley Publishing Co., Inc., Redwood 
City, California, ISBN 0-201-1 1603-0 

Krauss L.M., 1998, ApJ 501, 461 

Masuda N., Eriguchi Y., 1997, ApJ, 489, 804 

Masuda N., Nishida S., Eriguchi Y., 1998, MNRAS, 297, 1139 

Meszaros R, 2002, ARA&A, 40, 137 

Perlmutter S. et al. 1999, ApJ, 517, 565 

Rezzolla, L., Yoshida S'i., Maccarone T. J., Zanotti, O., 2003a, 

MNRAS, 344, L37 
Rezzolla, L., Yoshida S'i., Zanotti, O., 2003b, MNRAS, 344, 978 
Riess A. G. et al. 1998, Astron. J., 116, 1009 
Stuchlik Z., Hledik S., 1999, Phys. Rev. D, 60, 044006 
Stuchlik Z., Slany R, Hledik S., 2000, A&A, 363, 425 
Zanotti O., Rezzolla L., Font J. A., 2003, MNRAS, 341, 832 



